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ABSTRACT 

Star- forming galaxies at redshifts z > 6 are likely responsible for the reionization of the universe, 
and it is important to study the nature of these galaxies. We present three candidates for z ~ 7 
Lyman— break galaxies (LBGs) from a 155 arcmin^ area in the CANDELS/COSMOS field imaged by 
the deep FourStar Galaxy Evolution (zFourGE ) survey. The FourStar medium-band filters provide 
the equivalent of R ^ 10 spectroscopy, which cleanly distinguishes between z ^ 7 LBGs and brown 
dwarf stars. The distinction between stars and galaxies based on an object's angular size can become 
unreliable even when using HST imaging; there exists at least one very compact z ^ 7 candidate 
(FWHM~0.5-1 kpc) that is indistinguishable from a point source. The medium-band filters pro- 
vide narrower redshift distributions compared with broad-band-derived redshifts. The UV luminosity 
function derived using the three z ^ 7 candidates is consistent with previous studies, suggesting an 
evolution at the bright end {Muy ~ —21.6 mag) from z ^ 7 to z ^ 5. Fitting the galaxies' spectral 
energy distributions, we predict Lyman-a equivalent widths for the two brightest LBGs, and find 
that the presence of a Lyman-a line affects the medium-band fiux thereby changing the constraints 
on stellar masses and UV spectral slopes. This illustrates the limitations of deriving LBG properties 
using only broad-band photometry. The derived specific star- formation rates for the bright LBGs are 
^ 13 Gyr~^, slightly higher than the lower-luminosity LBGs, implying that the star-formation rate 
increases with stellar mass for these galaxies. 

Subject headings: galaxies: high-redshift — galaxies: Lyman break galaxies-galaxies: Luminosity Func- 
tion 



1. INTRODUCTION 

Discovering high-redshift {z > 6) galaxies, and un- 
derstanding their physical nature is of great importance 
in studying the early universe because these objects 
are likely sources f or reionization of the intergalactic 
medi um flGM: e.g. PlVenti et al.l lmol: iSalvaterra et al I 
12011'; 'Finkelstein et alj |2012D "at ^ > 7 (jFan et al.l 120021 : 
[Komatsu et al. 2011). With the advent of large near- 
infrared (NIR) detectors, significant progress has been 
made in discove ring g a laxies at z > 7 (e. g. Ive et al. 20061 
IVanzella et all 120111 : lOno et al.n2012l : iSchenker et al.l 
12012"), as the rest-frame UV light from these galaxies 
shifts to A >1 fiiR. 

There are two main techniques for finding high-redshift 
(^ ^ '^) galaxies : (1) the Lyman-break or dropout 
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meth o d (e.g. iSteidel et all 11995' '1999; Adel berger et all 
20041: iDickinson et all 12004: .Giavalisco et al.l 12004 : 



Ouchi et all 120041 : iBouwens et al.n2007l : iMcLure et all 
2010f ). selected based on strong absorption in their spec- 



tral energy distribution (SED) at wavelengths shortward 
of redshifted Lyman-a due to IGM HI absorption, and 
(2) the narrow-band imaging technique, which identifies 
galaxies with a strong, redshifted Lyman-a emission line 
using narrow- band fil ters (e.g. iM alhotra & Rhoads"2 002l: 
Rhoadsetall 120051: Hve et al.l 2006; Kashikawa etj^ 



2006; Ouchi et al. '2009bl: iHibon et al.l 120101 : ITilvi et~all 
2010 : Krug et al. 2012). A unique advantage of narrow- 
band selected galaxies is that these galaxies are more 
likely to be confirmed via spectroscopic followup be- 
cause they are pre-selected based on the strong Lyman- 
a emission line. On the other hand, LBGs selected 
via the broad-band dropout method are unbiased in 
their selection (unlike Lyman-a selected galaxies), and 
LBG surveys probe larger cosmological volumes due 
to broader filte rs. Currently there are many candi- 
date LBGs (e.g. 'Bouwens e t"al]l2007l : lOuchi et al.l 120091 : 
iMcLure et al. 2010) at redsh ifts much greater than z ^ 7 
(e.g. Bouwens et al. 2011; Bradl ev et al.! '2012':'Ya n et al.l 
2011; Finkelstein et al. 2012 : Oesch et al. 2012 ). 

With growing samples of 2; ^ 7 galaxies, we are able 
to calculate the rest-frame UV luminosity function (LF), 
which is a simple yet powerful probe of the early uni- 
verse. The luminosity function allows us to constrain the 
observed number density of star-forming galaxies from 
which we can estimate (with some assumptions) the num- 
ber of hydrogen- ionizing photons available to reionize the 
IGM. In addition, the observed luminosity function al- 
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Fig. 1. — Filter transmission curves of the FourStar medium- 
band compared to other broad-band filters. The top panel shows 
broad-band Subaru z\ HST ^los, -^125 ^ Hiqq filters, while the 
bottom panel shows FourStar medium-band filters Ji, J2, J3, Hs 
and Hi, which provide R~10 resolution spectroscopy. Black line 
shows a template spectral energy distribution of a 2; = 7.4 galaxy 
with the Lyman-o; emission line falling in the Ji filter. The top 
axis marks the Lyman-o; redshifts as a function of wavelength. 



lows US to study the galaxy evolution by comparing the 
number density of galaxies at different redshifts. 

Several studies have focused on understanding the evo- 
lution of the UV luminosity function from z ~ 3 to 6 
and find that there is a significant decrease in the num- 
ber density of UV bri ght galaxies from low to h igh red- 
shifts (e.g . Stanwav et al. 2003|; IShimasaku et al.l I2QQ5I : 
iBouwens et a l. 2006). If this trend continues, we ex- 
pect even a greater decline in the number density of 
galaxies at redshift z > 6. Current observations at 
z > 6 suggest a strong evolutio n in the cha racteris- 
tic m agnitude M"*" ^ -19.8 mag ("Bouwe ns et al.l 120081: 
McLure et al.ll2Q10l: lOuchi et al.l [2 009: Ya n et al.l IMol : 
Castellano et al.l 120101 : lBouwens~e t al. 201 Q|) compare d 



with M* --20.9 mag at z - TTBou wens et al.ll2007[ ). 
Most of these studies, especially at z ~ 7, rely on 
HST observations which have the advantage of prob- 
ing fainter galaxi es, however w ith smaller surve y areas . 
Recentlv. iQuchi et al.l (J2009D. ICastellano et al] (j2010f ). 
iCapak et all (j2011f ). and lBowler et al.l (j2012r ) found sev- 
eral candidates using ground-based observations, which 
have the advantage of surveying larger areas and thus 
probing brighter, rarer galaxies. 

While it is relatively easy to construct the luminosity 
function, its accuracy depends on a number of systemat- 
ics, including the degree of contamination of the z ^ 7 
galaxy sample. It is affected by contaminants that mimic 
the colors of high-redshift galaxies. The main contam- 
inants of z ^ 7 LBGs are Galactic brown dwarfs, low- 
redshift {z ^ 1 — 2) dust obscured galaxies with very 



faint continuum in the visible bands, and transient ob- 
jects. The brown dwarf stars may be separated from 
(resolved) galaxies using high angular-resolution imag- 
ing. However, the ability to resolve galaxies can become 
unreliable at the faint magnitudes and compact sizes typ- 
ical of z ^ 7 galaxies, even with the HST observations. 
For example, the rejection of point sources from z ^ 7 
samples precludes the possibility that some z ^ 7 galax- 
ies have morphologies that are unresolved even at HST 
resolution (FWHM ^0'.'2 which corresponds to 1 kpc at 
z = 7), leading to sample incompleteness. These uncer- 
tainties need to be accounted for when constructing the 
luminosity function. 

In addition to studying the UV luminosity func- 
tion, understanding the physical properties of high- 
redshift galaxies is important to link their evolu- 
tion during the early stages of galaxy formation 
to galaxies at other epochs. As we probe red- 
shifts z > 2, galaxies on average have smaller 
stellar masses ~ 10^^ M^t^ ( e.g. ISawicki fc Ye^ 
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Yan et al." ■2005'; "Erb et al. ■2006'; "Fontana et alj 
Reddv et al. 2006; Overzier et a l. 2009), bluer UVcol- 
ors (e.g. Dickinso n et al.l 120031: iPapovich et al.l 120041: 



Labbe et al.. .2007: lOverzier et al.l 120091 : IBouwens et all 
2009al) . and younge r stellar populations with smaller 



amounts of dus t (e.g. iMeurer et al."1999': 'Papovich et al.l 
I2QQTI: Shaplev 'etaLl I2001L 12005: Reddv et al. 20^ 
2006, 2008; Reddy & Steidel 2009; B ouwens e t al. 2009a; 
Finkelstein et al. 2010 . 2011) . Moreover, the specific 
star-formation rates (sSFR: the SFR per unit stellar 
mass) of galaxies at 2; > 2 are roughly constant as a func- 
tion of redshift, with sSFR ranging fro m ^ 2 to 6 Gyr~^, 
albeit with large u ncertainties (e.g. Stark et al. 2009; 
Gonzalez et al.l 12010 : Bouw ens et al . 2011; Reddy et al.: 

>T2l ). In all cases, the physical properties of the 
most distant galaxies {z > 6), are derived by fit- 
ting the spectral energy distribution templates to the 
galaxy photometry measured in broad-band filters (e.g. 
[Fmkelstein et al. 2011). While some effects such as the 
degeneracy among different derived parameters are al- 
ready known, the template fitting to the broad-band pho- 
tometry may introduce (as of yet) unknown systematic 
errors on the constraints of these parameters. 

The specific goals of this paper are (1) to use near- 
IR medium-band filters, for the first time, to reliably 
identify z ^ 7 galaxies and to compare the effective- 
ness of the medium-band technique to that of broad- 
band photometry; (2) to discuss the contamination rate 
of brown dwarfs in z ^ 7 galaxy samples; and (3) to 
derive the physical properties of the galaxies from the 
medium-band photometry and to test for systematic un- 
certainties on physical properties using broad-band data 
alone. Photometric redshifts derived using medium-band 
filters have much narrower probability distributions com- 
pared to the broad-band derived photometric redshifts 
owing to the higher spectral resolution (R~10) of the 
medium-band filters. This led to a recent discovery of 
a distant (z = 2.1) galaxy cluster in the COSMOS field 
(|Spitler et al.ll2012l ). Owing to very deep imaging and 
excellent seeing, we are also able to place better con- 
straints on the redshifts and on the evolution of bright 
star-forming galaxies out to z ^ 7. 
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Fig. 2. — Left panel: Curves-of-growth for point sources in the convolved and unconvolved images. Dashed and solid lines show flux 
ratio before and after convolution, respectively, where / is the flux in a given band and /^z is the z' (reference image) flux enclosed 
within radius r. Different colors indicate different fllters. The flux in aperture radii larger than 0(^48 among convolved images contribute 
to < 2% uncertainty for point sources. Right panel: The signal-to-noise (S/N) for point sources as a function of radii, in the convolved 
images. Different lines show the S/N of same stars (here shown only the J3 fllter for illustration purposes) that were used for creating the 
convolution kernel (Section 2). For photometry, we use a circular aperture with 0(^55 radius as this yields maximum S/N. 



TABLE 1 

Photometry for COSMOS field. 



Filter 


\{lim) 


Depth (mag^ 


Filter width 

(A)_ 


9' 


0.44 


27.96 


1384 


r' 


0.62 


27.95 


1460 


i 


0.83 


26.59 


2270 


z' 


0.90 


26.43 


1384 


Ji 


1.05 


26.10 


1030 


J2 


1.14 


26.07 


1410 


J3 


1.28 


25.69 


1320 


Ji25(F125W) 


1.25 


26.93 


2991 


Hs 


1.55 


25.15 


1600 


Hi 


1.70 


25.17 


1610 


iJi6o(F160W) 


1.54 


27.04 


2880 


K, 


2.15 


25.13 


3380 


[3.61^ 


3.56 


26.0 


7274^ 


[4.5^ 


4.52 


26.0 


9914^ 



^ 5cr depth in 1^1 diameter aperture 

^ 3cr depth in 2^4 diameter aperture (Ashby et al. 2013 sub- 
mitted) 
^ Filter width where fllter transmission is >50%. 
^ Transmission is >30%. 



This paper is organized as follows. Section 2 describes 
the observations, data reduction, and catalog generation. 
In Section 3 we present our LBG color selection method 
and sources of contamination and their identification. 
Section 4 describes the construction of the UV luminos- 
ity function and its comparison to previous studies. We 
compare the advantage of using medium-band over only 
broad-band photometry, and its influence on the derived 
physical properties of high-redshift galaxies in Section 5. 
We summarize our conclusions in Section 6. Throughout 
this paper we use AB magnitudes and standard cosmol- 
ogy with h = 0.7, fix = 0.7, and Qm = 0-3. 

2. OBSERVATIONS AND DATA REDUCTION 



We observed a 155 arcmin^ area in t he CAN- 
DELS/COSMOS field (jScoville et al.l 120071) centered 
at RA: 10:00:31 and Dec:+02:17:21, using the FourStar 
instrument, as part of the ongoing FourStar Galaxy 
Evolution survey (zFourGE during 2011-2012. The 
FourStar instrument (Persson et al. 2013) is a near-IR 
camera mounted on the Magellan/Baade telescope. It 
has 4096 x 4096 pixels with a 0''.'159/pixel resolution 
and covers a 10^8xl0^8 field of view. We observed 
with five adjacent medium-band filters (Ji, J2, J3, Hg, 
Hi'. Table 1) with a resolution of i? ~ 10, and one 
broad-band {Ks) filter (Labbe et al in preparation). Fig- 
ure 1 shows filter transmission curves for the FourStar 
medium-band filters compared to the broad-band filters 
from Subaru/Suprime (z'), and HST (F105, -^125, ^i6o) 
filters. The Ji filter has a similar central wavelength 
to the Fios filter (and also to many F-filters on ground- 
based telescopes), but the Ji filter is much narrower com- 
pared with the HST F105 filter. 

We processed all raw images using a custom de- 
signed IDL-based pipeline (described in more detail in 
Labbe et al in preparation) adapted from the NEW- 
FIRM Medium Band S urvey data reduction pipeline 
(NMBS: Whita ker et"aIll2Qllf ). The data quality of the 
FourStar images is excellent and the point-spread func- 
tion (PSF) FWHM for the stacked FourStar images cor- 
respond to 0.55, 0.52, 0.51, 0.51, 0.57, & 0.48 arcsec in 
the Ji through Ks filters, respectively. In addition to 
this data set, we used publicly available data in the op- 
tical (Subaru ^', r', z' ^ & ACS ^8i4(^))5 deep near-IR 
HST WFC3 Ji25 fc H^m from the CAN DELS survey 
(iKoekemoer et a02011l:lGrogin et al.ll201lt ) , and infrared 
IRAC ("FazioelalJ |200J) data. We also used deep 3.6 
and 4.5/im images (Ashby et al 2013 submitted) from the 
Spitzer Extended Deep Survey {SEDSf^, 5.8 and 8.0/im 

^^ http://candels.ucQlick.or g/ 

^^ |http://z-fourge. strw.leidenuniv.nl/" 
http 
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dat43 from the Spitzer Space Telescope (| Werner et al.l 
|2004) . Table 1 presents the central wavelengths and lim- 
iting magnitudes for the zFourGE data as well as infor- 
mation for the ancillary data. 

The PSF among different bands varies due to differing 
central wavelengths, different instruments used for imag- 
ing, and in the case of ground-based observations, due to 
different seeing conditions. We convolved all images to 
a reference image having the broadest PSF-FWHM (in 
our case the Subaru z' image with PSF FWHM=0y8). 
In each band, we first constructed an empirical PSF by 
identifying stars and stacking them with their centers 
aligned. We then used an I DL routine, DECONV_TOOL 
(jVarosi fc Landsmanll 19931 ). to construct a kernel to con- 
volve the image PSF to match the stellar PSF. Figure 2 
shows the curve-of-growth of point sources in the images 
following the convolution process. As shown in the fig- 
ure (left panel), the relative fiux for point sources in each 
band is matched to better than 2% for all apertures with 
radii larger than 0V47. 

The photometry in circular apertures was optimized 
for faint, compact sources. Figure 2 (right panel) shows 
the S/N for point sources measured in circular apertures 
as a function of aperture radius. The S/N is maximized 
for point sources with apertures of radii 0^55. Therefore 
we use circular apertures with diameter of V.\ for the 
photometry in our analysis. 

2.1. Spitzer/IRAC Photometry 

As mentioned earlier, to extend the wavelength cover- 
age of our candidates and to derive the physical proper- 
ties as accurately as possible, we used recent deep Spitzer 
IRAC data (Ashby et al 2013 submitted) in two bands 
(3.6 /iTTi, 4.5 jiim) obtained from the SEDS survey. How- 
ever, due to its low resolution, neighboring objects tend 
to blend resulting in incorrect fiux mea surements. To cir- 
cumvent this issue, we used GALFIT (Peng et al.n2002') 
code to measure the object fiuxes, as described below. 
First, we fit the two dimensional surface brightness pro- 
files of all objects within 7^' of the z ^ 7 candidate (ex- 
cluding the z ^ 7 galaxy) and subtracted these fits to 
obtain a residual image. We then measured photometry 
of the z ^ 7 candidates from the IRAC data using the 
IDL procedure APER, measuring the fiux in a fixed 2" 
diameter aperture and correcting to the total fiux using 
the aperture corrections derived by taking the ratio be- 
tween 2^' and lb'' diameter aperture fiuxes of isolated 
sources. Finally, we derived the uncertainties by mea- 
suring the rms in apertures of the same size as above, 
randomly placed in regions devoid of objects in the IRAC 
images. 

2.2. Source Catalogs 

We identified sources in each of the ima ges using Source 
Extractor software (SExtr actor : Bertin fc ArnoutsI 
119961 ) in dual image mode. In this mode, a detection 
image is used to identify the pixels associated with each 
object, while the fiuxes are measured from a distinct pho- 
tometry image. For the detection image, we used a x^ 
image, which optimally sums multiple ima ges (account- 
ing for varying S/N) for source detection (|Szalav et al.l 

1^ |http://irsa.ipac.caltech.edu/data/COSMOS7] 



|1999[ ). The x^ image {Jchq) w as constructed from the 
Ji, J2, J3 images using SWarp (jBertin et al.ll2QQ2[ ). We 
measured object photometry from the PSF- matched im- 
ages in a fixed aperture size of 1.1^' diameter. 

While SExtractor estimates the noise for each object, 
this estimate is incorrect in the convolved images as it 
does not account for correlated noise, which results from 
the image reduction and convolution processes. To esti- 
mate the fiux uncertainty for each object in the convolved 
image s we fo llowed an empirical approach described in 
Labbe et al.l ([2b03). To estimate the pixel-to-pixel noise 
we measured the rms background variation as a func- 
tion of linear size N = a/A where A is the area of an 
aperture. First, we measured the sky fiux in increasingly 
larger apertures that are randomly placed in the PSF- 
matched-convovled image, taking care to avoid detected 
objects. From the distribution of fiuxes measured in a 
given aperture size, we measured the standard deviation, 
<j. The 1(7 erro r on the flux was then calculated using 
equation 5 from IWhitaker et al.l ()2Qllf ). Our final cata- 
log contains aperture fluxes obtained from PSF-matched 
images in g' , r' , i, z' , Ji, J2, J3, J125, Hs, Hi, i^ieo, and 
Ks along with fiux uncertainties. 

2.3. Photometric Redshifts 

We derive photometric redshifts using the ph otomet- 
ric redshift code EAZY (Brammer et ajj l2QQ8f ) taking 
advantage of extensive multi-wavelength data publicly 
available in the COSMOS field for finding the best-fit 
solution to the obser ved data. The EAZY code (v2.1; 
iBrammer et al.ir2Qllf ) contains seven different templates 
spanning a large range of galaxy colors and incorporates 
nebular emission lines. While our z ^ 7 candidate se- 
lection is based primarily on a color-color selection, to 
increase the reliability of our candidates, we also require 
that each of the candidates have photometric redshifts 
z>6. 

3. LBG COLOR SELECTION AT Z ^ 7 

Our initial color-selection criter ia to identify candid ate 
galaxies at z ~ 7 is adapted from lOesch et al.1 (|2Q1C)[ ): 

SmJchq) > 7.0 (1) 

S/N(optical) < 2.0, 

Ji — J3 < 1-3 mag, 

z — Ji > 0.8 mag, 

z - Ji> 0.9 + 0.75( Ji - J3) mag, 

z' - Ji> -1.1 + 4( Ji - J3) mag. 

Similar criteria have been used in other studies to iden tify 
z ^ 7 galaxies (e.g. Bouwens et al. 2009a; Ouchi et al.l 
120091 : ICastellano et al. 2010). In addition to these crite- 
ria, we required that each of the candidates be detected 
in the J2 filter as this filter is very effective in discriminat- 
ing T-dwarf stars (Tilvi et al in preparation ) , due to their 
methane and water absorption features (jTinnev et al.l 
'2012), from high-redshift star- forming galaxies. 

Figure 3 shows a J1-J3 vs z'-Ji color-color plot and 
illustrates the region (dashed line) that corresponds to 
the selection criteria of Equation 1. Color-color tracks 
of star-forming galaxies with two different dust attenu- 
ation are shown in blue and red lines with redshift la- 
bels. These tracks are obtained from .Bruzual fc CharlotI 
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Fig. 3. — Color-color plot for z ~ 7 LBG candidate selection. Filled circles, triangle and square symbols indicate all the objects that 
pass the LBG selection criteria (Equation 1; dashed line), adapted from Oesch et al. (2010). Blue and red lines show the expected colors 
of a star-forming stellar population, for two extinction values, as a function of redshift (redshifts indicated for reference) using models 
from Bruzual & Chariot (2003). Star symbols show the expected colors of very late-type stars derived from observed spectra of M, L, 
and T-dwarf stars. The initial color selection (dashed line) yields four z ~ 7 LBG candidates (filled circles), four brown dwarf stars (filled 
squares) and one likely foreground galaxy (filled triangle). Based on the galaxy tracks (blue and red lines ) and brown dwarf star colors, 
we refined the z ~ 7 LBG selection criteria (Equation 4; shaded region) that eliminates most of the nearby stars. This new selection yields 
three robust z ~ 7 LBG candidates. 



()2QQ3f ) stellar population synthesis models by integrating 
the model spectral energy distributions with the z' ^ Ji, 
and J3 filter. The models assumed E(B-V)= 0.0 (blue 
line) and 0.1 (red line) with Z=O.2Z0 metallicity, con- 
stant sta r-formation history, IGM opacity taken from 
(|Meiksinl i2006). and a constant age of t= 100 Myr. Based 
on these color tracks it can be seen that the star-forming 
galaxies occupy a specific region of the color-color space. 
We identified 9 objects in total satisfying the selec- 
tion criteria from Equation 1. Based on several tests de- 
scribed in the following sections, out of the 9 objects in 
the initial sample (hereafter initial source list), four ob- 
jects are identified as z ~ 7 LBG candidates, four objects 
are identified as brown dwarf stars and the remaining 
object is identified as a likely foreground dust-obscured 
galaxy at z < 3. To minimize the contamination from 
stars and foreground galaxies, we further refine the color- 
color selection (Section 3.3) that include only the robust 
z ^ 7 LBG candidates (Figure 3; shaded region). Figure 
4 shows postage stamps of these z ^ 7 LBG candidates, 
as well as one likely T-dwarf star for comparison. 

3.1. Sources of Contamination 

The color-color selection of z ^ 7 LBG candidates suf- 
fers from several possible sources of contamination in- 
cluding image artifacts (or spurious detections), nearby 
cool brown dwarf stars which mimic high-redshift galaxy 
colors, intermediate redshift galaxies with faint contin- 



uum in the visible bands, and transient sources including 
super novae or asteroids. In this and the following sec- 
tions we discuss in detail the likelihood of any of these 
contaminating sources entering the initial source list. 

3.1.1. Spurious detections 

In the z ^ 7 LBG selection, we require that each of the 
objects is detected with S/N > 4 in more than one band 
redward of the z^ filter. This requirement makes spurious 
detections due to either background or electronic noise 
extremely unlikely. 



3.1.2. Transient/ High-proper motion objects 

It is possible that transient objects including asteroids, 
and supernovae can be detected in more than one band. 
However, with the availability of two epoch data, sep- 
arated by about one year, we can eliminate any tran- 
sients by comparing their positional shift between the 
two epochs. Therefore, we conclude that the initial 
source list does not contain any transient sources. 

3.1.3. Foreground galaxies 

While the z ^ 7 candidates are required to have no 
detection (< 2<j) in the observed visible bands, it is 
possible that the initial source list contains foreground 
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Fig. 4. — Image cutouts of four z r^ 7 LBGs and one of the T-dwarf candidates (zF26680). Top panel shows color composite {z' ,Ji,Hiqo) 
image cutouts while lower two panels show image cutouts in the individual filters. The four z ~ 7 LBG candidates (zF13965, zF25025, 
zF26836, zF26343) ) are selected from the color-color selection (Equation 1) while zF26343 is excluded in the refined color-color selection 
(Equation 4). Each candidate is marked within a 2'' radius circle. These candidates are clearly detected in at least one of the medium-band 
filters while undetected (< 2a) in any of the visible bands. The integrated probability distribution Vq (Equation 2) is shown on the 
right side. In addition to satisfying the color-color selection criteria, each of the candidates is required to have 7^6 > 0.7 and S/N< 2 in 
the visible x^ image. These two criteria minimizes the contamination from low-redshift galaxies. In the bottom panel, we show one of 



the candidate T-dwarf stars that was selected in the initial color-color selection (Oesch et al 2010). 
distinguished from galaxies using medium-band photometry. 



Such T-dwarf stars can be cleanly 



objects, likely z ^ 1 — 3 highly dust obscured galaxies 
that have very low fluxes in the visible bands, but are 
bright in the near-IR and IRAC bands, (see Labbe et al 
2005, Papovich et al 2006). These low-redshift galaxies 
cannot be completely ruled out, but the contamination 
from such objects can be minimized because they oc- 
cupy a specific color-color space (Figure 3) compared to 
the z ^ 7 galaxies which have bluer colors. Further- 
more, many of the red, intermediate redshift galaxies at 
z ~ 1 — 4 have bright IRA C fiuxes with ty pical colors 
[4.5] - [5.8] < AB mag fe.g. lYan et all[200l . which dis- 
tinguishes them strongly from z ^ 7 galaxy candidates. 
In addition to these signatures, the photometric redshifts 
can provide an additional constraint on the likelihood of 
an object being a low-redshift galaxy. 

In order to compute the likelihood of our candidates 
being z ^ 7 galaxies, we followed a similar procedure, as 
described in Finkelstein et al. (2010). We compute the 
total redshift probability distribution Vq given by 



V6 = 



IeP(z)dz 
J^P{z)dz' 



(2) 



where Jg P{z)dz is the integrated probability distribu- 
tion from z = 6 to z = 9. This probability distribution 

is normalized by J^ P{z)dz such that Vq < l- We make 
use of the probability distribution P{z) derived using the 



photometric redshift code EAZY and require that each 
of the z ^ 7 candidates has Vq > 0.7. This means that 
each of the z ^ 7 candidates has greater than 70% prob- 
ability of being a z > 6 galaxy, based on the photometric 
redshift. All of the 9 objects in the initial source list 
have Vq > 0.99 and in principle, we could restrict the 
sample to the stricter requirement of Pe > 0.99, but this 
would have risked missing some z ^ 7 gala xies near our 
detection limit (see iFinkelstein et al.ll201Qf ). 

3.1.4. Visible ^ Image 

To further test the reliability of z ~ 7 candidates' non- 
detection in the individual visible bands, we combine g' , 
r', and i images to construct a visible x^ image; the 
combined x^ image is significantly deeper than individ- 
ual images. We then run SExtractor to measure the S/N 
of all objects in the x^ image. Based on S/N(> 2) in 
the visible x^ image, we reject one object (zF4329) out 
of 9 objects from the initial source list, as being a likely 
foreground galaxy with very faint continuum. This re- 
jection i^_also_su2Ported by its strong detection in the 
MIPS (|Rieke et al.ll2004f ) 24/im fiux; there are no MIPS 
detections for any of the remaining objects. 

3.2. Distinguishing Stars from Compact Galaxies 

Galactic brown dwarfs (M, L, and T ) are one of the 
main contaminants because their near-IR colors resem- 
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ble high-redshift {z ^ 7) galaxy colors. While it is pos- 
sible to identify candidates for stars using their FWHM 
and stellarity index (produced by SExtractor) when us- 
ing HST observations, these classifiers tend to become 
unreliable at fainter magnitudes. In addition, at high- 
redshifts we expect some galaxies to be very compact 
with their FWHM (< 1 kpc at 2: ~ 7) comparable to 
the PSF FWHM. In such cases, even with HST obser- 
vations, it would be challenging to discriminate between 
stars and compact galaxies. 
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Fig. 5. — FWHM as a function of magnitude. Stars (star 
symbols) are identified based on their FWHM and stellarity 
index >0.8 extracted from SExtractor. Filled symbols (same 
as in Figure 3) show all objects that passed the initial 2; ~ 7 
color selection. Top panel: J2 mag vs FWHM from the FourStar 
bands (ground-based observations). Bottom panel: J125 mag vs 
FWHM (space-based observations). The distinction between stars 
and compact high-redshift galaxies, based on their FWHM, is 
unreliable from ground-based observations, compared with HST 
observations. All fours dwarf stars (square symbols) have stellarity 
indices > 0.7 while object zF25035 has a stellarity index of 0.6 
despite being its FWHM similar to the stellar FWHM. Small dots 
represent all other objects in the catalog. For clarity purposes, 
we show only a few representative objects. Note that the J125 
magnitude for the foreground object (triangle) is significantly 
fainter compared with the J2 magnitude due to the fact that 
in the raw (0(^06) J125 image, this object is identified as three 
different objects. 



In the following subsections, we perform three tests to 
distinguish stars from compact galaxies: (i) tests com- 
paring the objects' FWHM and stellarity indices, (ii) 
tests of the objects' spectral features indicative of molec- 



ular absorption in stellar atmospheres of brown dwarf 
stars, and (iii) tests of the objects' surface brightness 
profiles. 

3.2.1. Identifying Stars using FWHM and Stellarity Index 

Figure 5 shows the FWHM (obtained from SExtractor) 
for objects as a function of J2 magnitude as well as the 
FWHM measured in the (unconvolved) HST J125 band 
to compare ground-based data with the space based ob- 
servations. The horizontal dashed line in each panel of 
Figure 5 indicates the stellar PSF FWHM obtained by 
stacking ten isolated stars in the respective images. 

In addition to the FWHM, SExtractor produces a stel- 
larity index (CLASS STAR), which determines a likeli- 
hood of an object being either a point or an extended 
source using the self-training of a neural network. A stel- 
larity index closer to 1 is likely a point source while ex- 
tended objects like galaxies have stellarity indices closer 
to 0. Figure 5 (top panel), shows some objects classified 
as stars based on their high (> 0.8) stellarity index (star 
symbol). All objects in the initial source list are shown 
in filled circle, box, or triangle symbols. As can be seen 
from the top panel, most of the objects appear to be re- 
solved in the FourStar image. However, this appears to 
be due to the seeing conditions and the magnitude limit 
of the ground-based data. The lower panel in Figure 5 
shows that five of the objects in our sample (zF2014, 
ZF24834, ZF30058, zF26680, and zF25035) have FWHM 
in the HST image much closer to that of a point source. 
Four of these objects (square symbols in Figure 5) have 
high stellarity indices (>0.7), indicating that they are 
likely point sources. In fact, fitting brown dwarf spec- 
tral templates to the medium-band observations (Section 
3.2.2), we conclude that all four of these are candidates 
for brown dwarf stars. The remaining compact source, 
zF25035 has a stellarity index of 0.6, suggesting that it 
might be a compact galaxy rather than a point source. 
In the following sections, we discuss additional tests to 
distinguish between the possibility that this object is a 
point source, or a very compact z ^ 7 galaxy candidate. 

3.2.2. Spectral Fitting to Medium-band Photometry 

Due to the availability of medium-band photometry, in 
addition to using stellarity indices and FWHM, we can 
use a spectral fitting technique whereby observed spec- 
tral templates are fitted to the medium-band photome- 
try. This method works very well especially in identify- 
ing T-dwarf stars because the central wavelengths of the 
FourStar medium-band filters trace the strong absorp- 
tion features of T-dwarf stars; this would not be pos- 
sible u sin^ broad-band photometry. Ivan Dokkum et aP 
(|2009l ) have shown one of the ways to identify T-dwarf 
stars using medium-band photometry from NEWFIRM 
(see their Figure 5). 

To identify dwarf stars, we synthesized medium-band 
pho tometry from the M, L , and T-dwarf observed spec- 
tra (jBurgasser et al.l [20061 ) and fit these to the FourStar 
photometry of all candidates to obtain the minimum x^ 
between the candidate photometry and the dwarf tem- 
plate spectra. Figure 6 demonstrates the effectiveness 
of medium-band filters in identifying T-dwarf stars and 
best-fit spectral energy distribution (SED) of z ^ 7 
LBGs. The bottom-left panel shows the best-fit spec- 
trum of an observed T-dwarf star, which clearly traces 
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Fig. 6. — Spectral energy distributions (SEDs) of the three z ^ 7 LBG candidates and one of the T-dwarf candidates. For all panels the 
filled circles represent observed photometry. For the three panels showing the 2; ~ 7 LBG candidates (zF25035, zF13965, zF26836) the 
solid line shows the best-fit SED template and the inset panel shows the photometric probability distribution obtained from the EAZY 
code. In the left panels, also shown are the image cutouts of the LBG candidate (top left) and T-dwarf candidate (lower left panel). 
For one of the the T-dwarf candidates (lower left panel), the solid line represent the best-fit brown dwarf template spectrum, derived by 
minimizing x^ between the Four Star photometry and the expected fiux density synthesized from spectral templates of M, L, and T-dwarf 
stars from Burgasser et al. (2006). The medium-band photometry clearly traces the T-dwarf observed template spectrum (lower left 
panel) and can distinguish between different T-dwarf types (Tilvi et al., in preparation). 



the medium-band photometry, while the remaining pan- 
els show z ^ 7 candidates fitted with galaxy SEDs. Using 
this method, we classify four objects, previously identi- 
fied as stars based on their J125 FWHM, as dwarf stars. 
For the compact object (zF25035) however, we did not 
find a good fit with any of the M, L, or T-dwarf observed 
spectra. Thus, this test favors the conclusion that this 
object is a compact z ^ 7 galaxy. 

3.2.3. Surface Brightness Profile 

Figure 7 shows the surface brightness (J125) profiles 
of five objects: two T-dwarf candidates (with medium- 
band colors indicative of brown dwarfs and high stellar- 
ity); two z ^ 7 candidate galaxies including the com- 
pact zF25035, and a random object selected from the 
subsample of all objects with stellar ity > 0.9. The two 
T-dwarf stars both have stellarity indices >0.9 and have 
FWHM nearly equivalent to the expected WFC3 PSF 
FWHM (OV2). Their surface brightness profiles are also 
indistinguishable from each other. While the compact 
z ^ 7 galaxy has a FWHM equal to that of the point 
sources, its surface brightness profile resembles the other 
z ^ 7 candidate galaxy. This is also supported by its 
lower stellarity index of 0.6 and thus strengthens our 
earlier conclusion that this object (zF25035) is a com- 
pact z ^ 7 LBG candidate. Moreover, it shows that at 



2: ~ 7, there is a population of very compact (FWHM <1 
kpc) star-forming galaxies that are unresolved in typical 
image quality of ground-based imaging. Note however, 
that the FWHM or surface brightness profile can become 
unreliable at fainter magnitudes. 

To demonstrate this, we have chosen a random star-like 
object with stellarity > 0.9 and FWHM similar to those 
of the compact galaxy, but with slightly fainter magni- 
tude. As illustrated in Figure 7 (grey dashed line), the 
surface brightness profile of the random object is indis- 
tinguishable from galaxies despite its stellar-like FWHM. 
This comparison demonstrates that the star versus com- 
pact galaxy classification based on either FWHM, stel- 
larity index, or even surface brightness profile would 
lead to mis-identification, especially at fainter magni- 
tudes {H > 25 mag), even with HST-like image quality. 
In such cases, the near-IR medium-band photometry will 
allow us to cleanly distinguish between the spectral fea- 
tures of brown dwarf stars and z ^ 7 galaxy candidates. 

In summary, out of five compact objects, four objects 
are classified as brown dwarf stars, and the remaining 
object zF25035 is classified as a compact galaxy. Elim- 
inating such objects due to their misidentification could 
affect the bright end of the UV luminosity function where 
objects are rare. This may especially be of concern for 
the case of ground-based data when such compact galax- 
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ies are unresolved. 

3.3. Final Sample and New Selection Criteria using 
Medium-bands 

Using the initial color selection criteria defined in 
Equation 1, out of the 9 objects in the initial source list, 
four objects are identified as nearby dwarf stars based on 
their FWHM and spectral template fitting test, one ob- 
ject is likely a foreground galaxy based on the S/N(> 2) 
in the visible x^ image, and the remaining four objects 
are identified as z ~ 7 LBG candidates. To select the 
most robust z ^ 7 LBG candidates, we revise the z ^ 7 
candidate selection criteria using medium-band filters as 
follows. 



S/N{Jchg) > 7.0 (3) 

S/N(J2) > 4.0 

SmxLme) < 2.0, 

Ji — Js < OA mag, 

z' - Ji > 0.8 mag. 

This color-color selection is indicated in Figure 3 as the 
shaded region within the dotted polygon. As can be 
seen, this color-color selection isolates z ^ 7 galaxy can- 
didates from both low-redshift red galaxies and brown 
dwarf stars. Therefore, our new color selection criteria 
using the FourStar medium-band filters appears to pro- 
vide a z ~ 7 LBG sample that is free of contaminants, 
although future surveys (including spectroscopic confir- 
mations) are needed to confirm this. The photometry of 
all three z ^ 7 candidate galaxies is shown in Table 1, 
while Figure 4 and Figure 6 show postage stamps and 
best-fit SEDs obtained from EAZY, respectively. 

While this refined selection yields three robust z ^ 7 
LBG candidates, one of the candidates (zF26343 in Fig- 
ure 3) is excluded by this selection (which would have 
been selected if we had used the selection criteria from 
Equation 1), and will lead to sample incompleteness 
thereby affecting the UV luminosity function. We pre- 
vent this incompleteness by simulating the effective sur- 
vey volume that is enclosed within the shaded region in 
Figure 3. 

4. LUMINOSITY FUNCTION AT Z - 7 

The UV luminosity function provides a direct obser- 
vational measurement on the number density of star- 
forming galaxies at a given redshift. We derive the 
completeness-corrected luminosity function (j){m) as 



^(m) 



1 N{m) 
dm Feff 



(4) 



where N{m) is the number of detected objects with mag- 
nitude between m and m + dm^ and Veff is the effective 
survey volume. 

4.1. The Effective Survey Volume 

In order to accurately estimate the observed luminos- 
ity function we must calculate the effective survey vol- 
ume since this volume depends on both the apparent 
magnitude of the selected objects and their redshift dis- 
tributions. We estimate the effective volume (Kff) as 
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Fig. 7. — Comparison of surface-brightness profiles (J125) of the 
compact z ^ 7 galaxy candidate (zF25035) and point sources with 
stellarity index >0.9. Blue and red lines show surface brightness 
profiles of z ~ 7 LBG candidates and T-dwarf candidates, respec- 
tively. Dot-dashed line is a random object chosen to have similar 
magnitude and FWHM as those of the compact galaxy, except 
that this random object has stellarity index > 0.9 identifying it 
as a star. If we were to identify this random object based on its 
FWHM or stellarity index it would be classified as a star, even 
though its surface brightness profile resembles a non-stellar object. 
With typical ground-based image quality, and even with HST 
image quality, discrimination between galaxies and stars based on 
their stellarity index is unreliable especially at fainter magnitudes. 




Kff(^) 



-p{m^ z) dz, 



(5) 



where p(m, z) is a probability of detecting a galaxy of 
apparent magnitude m at redshift z specific to the image 
quality, depth, and passbands of our dataset, and dV/dz 
is the comoving volume in a redshift interval dz. 

In addition, due to photometric errors, objects can 
scatter in and out of the color-color selection region 
(shaded region; Figure 3). In the following sections we 
describe how we estimate the effective volume of our sur- 
vey. In short, we first compute p(m, z) by inserting and 
then recovering artificial galaxies in the real images as a 
function of apparent magnitude and redshift. We then 
calculate dV/dz for each redshift bin. 

4.1.1. Modeling the Colors of High- z Calaxies 

To generate artificial galaxies that are representative 
of high-z galaxies, we created a grid of spectral energy 
distributions (SED s ) of s tar-forming galaxies using the 
iBruzual fc Charlotl ()2003f ) models. In these models we 
used the Chabrier initial mass function, solar metallic- 
ity, and exponentially rising star-formation history. For 
a given redshi ft, we si n iulate d the effect of dust and 
IGM using the iCalzettil ()1997[ ) starburst model and the 
iMeiksinI (|2Q06f ) IGM prescription, respectively. We then 
convolved the resultant SEDs to calculate the bandpass 
magnitudes and colors through the z' and Ji filters at 
various redshifts (see Figure 3). 

To calculate the J1-J3 colors, we followed a different 
procedure instead of using galaxy tracks because the 
model colors are nearly constant (see Figure 3) over the 
redshift range we are currently probing. Here we as- 
sumed that the flux at a given wavelength can be de- 
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TABLE 2 




















Photometry of candidate galaxies at z 


-7. 










ID 


RA 


Dec 


z 


Ji 


Ji 


Js 


-^^125 


Us 


Hi 


-^160 


Ks 


[3.5] 


[4.5] 


ZF13965 


150.12579 


2.26662 


>27.4 


25.15 


25.24 


25.32 


25.13 


25.41 


25.25 


25.06 


24.87 


24.04 


24.80 








- 


±0.06 


±0.07 


±0.22 


± 0.01 


±0.56 


±0.50 


±0.01 


±0.42 


± 0.20 


±0.20 


ZF25035 


150.09906 


2.34362 


>27.4 


25.32 


25.39 


25.17 


25.33 


24.86 


25.47 


25.37 


25.24 


25.58 


23.64 








- 


±0.06 


±0.06 


±0.15 


±0.02 


±0.30 


±0.56 


±0.02 


±0.45 


± 0.50 


±0.20 


ZF26836 


150.08658 


2.31785 


>27.4 


26.39 


26.09 


26.28 


26.68 


27.09 


27.06 


26.47 


27.45 


>26.0 


>26.0 








- 


±0.14 


±0.11 


±0.41 


±0.05 


±2.26 


±2.31 


±0.11 


±3.41 


- 


- 


zF26343^ 


150.07767 


2.35306 


>27.4 


25.79 


25.72 


25.25 


25.35 


25.45 


24.91 


25.02 


25.11 


24.35 


24.42 








- 


±0.11 


±0.10 


±0.21 


±0.02 


±0.65 


±0.43 


±0.05 


±0.51 


±0.20 


±0.20 



All magnitudes are total magnitudes. 

Limiting magnitudes (2cr limits in I'.'l diameter aperture): 
^ This candidate falls outside the refined color selection. 



=28.9; r'=28.9; i=27.6. 



scribed as f\ ex A^ where /3 is the UV continuum slope. 
For a given Ji magnitude (calculated from the model 
SEDs as described in the previous paragraph), we calcu- 
late the J3 magnitude assuming j3 = —2.0. This value 
of [3 is consistent with z ^ 7 galaxies (JFinkelstein et all 

4T.2. The Probability Function p(m, z) 

We computed the probability function (completeness 
function) by inserting artificial sources based on their 
redshifts and colors, and recovering these artificial ob- 
jects using the same procedure as for the real galaxies. 
Now that we have magnitudes and colors of high-z galax- 
ies as a function of redshift, we first have to create an 
artificial object with some shape, and assign it a magni- 
tude. To create an artificial galaxy with its shape as close 
as a real observed galaxy, we choose a random object in 
a given magnitude bin from the science image itself (in 
this case Ji). 

To insert this object in the z' ^ Ji, and J3 images we 
choose 500 random positions avoiding any already exist- 
ing sources (> 3cr fiux) in the Ji image. Each object was 
assigned a z', Ji, and J3 magnitude based on the slope 
/3, and the color obtained from SED models, respectively. 
We simulate artificial sources in a redshift bin of 0.2 and 
magnitude bin dm = 0.1 mag. We then compute the re- 
covered fraction by taking the ratio of recovered sources 
(using the selection criteria as described in Section 3.3) 
to the inserted number of objects. This completeness 
accounts for both the detection incompleteness as well 
as the imperfect object selection due to photometric er- 
rors leading to additional incompleteness. To increase 
the number of simulated galaxies without increasing the 
crowding of sources in a single image, we repeat this pro- 
cedure ten times for a given redshift and magnitude bin 
and then average the recovered fraction. 

As an illustration we show p(m,z) for Ji=24.8 mag 
in Figure 8 (shaded area). To compare this complete- 
ness function to the broad-band filters, we also overplot 
the completeness function for the z y 7 selection using 
the ^850, ^105: and Jt25 fi lters from lOesch et al.l ()2012[) 
and iBouwens et al.l ()2012l ). We clearly see a broader 
completeness function for broad-band filters compared 
to our medium-band selection. In contrast, at Ji=24.8 
mag, p{m^ z) is much lower than that derived from much 
deeper HST broad-band surveys. Now that we have 
p(m, z), we calculate the effective volume by integrating 
the product of p(m,z) and dV/dz (equation 5). 

4.2. Evolution of UV Luminosity Function at z ^ 7 
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Fig. 8. — Probability of recovering artificial galaxies in our 
survey, shown for a source with Ji=24.8 mag (filled curve). 
Dashed and dotted lines show HST broad-band completeness 
function from Oesch et al (2009) and Bouwens et al (2012), 
respectively. The broad-band completeness distribution achieves a 
higher peak p{ni, z) due to the HST survey depths, but is slightly 
broader compared to the FourStar completeness distribution. The 
narrower distribution of medium-bands allows for much tighter 
redshift probability distribution. 



The apparent magnitude can be converted to the rest- 
frame UV magnitude Muv (at A = 1500 A) using 

M1500 = J3 - DM{z) - /c, (6) 

where DM is the distance modulus and k is the k- 
correction between the rest-frame wavelength 1500 A and 
J3 filter. Assuming a negligible k correction, the above 
equation can be rewritten as 

M1500 = J3 - 5 \og{DL/l pc) + 2.5 log(l + z) + 5, (7) 

where Dl is the luminosity distance at redshift z. 

We divided our z ^ 7 sample in two magnitude bins, 
with each bin 1 magnitude wide. The brighter luminos- 
ity bin contains two objects, while the fainter bin con- 
tains the remaining one object. Figure 9 shows the UV 
luminosity function of z ^ 7 galaxies derived from our 
observations for these three galaxies (blue filled circles). 

Currently there are only a few other studies explor- 
ing the UV luminosity f unction at magnitudes brighter 
than Muv= -21.5 mag (lOuchi et al.ll2009l: ICapak et all 
I20TTI : iBowler et al.l 120121 : iHathi et oim om. Our result 
is consistent with previous conclusions about the evolu- 
tion of the UV luminosity function: the UV luminosity 
function evolves from z ~ 5 to z ~ 7 at the bright end 



Discovery of Lyman Break Galaxies at 2; ~ 7 from the zFourGE Survey 



11 



-2 



o) -3 

05 



o -4 

Q. 



O 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

^ QzFourGE o zFourGE (+contaminants) 

♦ Ouchi+ 09 

A Bouwens+ 10a ^ " ^ 

L * Castellano+ 1 ^^5^ |^ ^ j- ■ i— y^ 



-5 

-6 

-7 



Bowler+ 12 
T McLure+ 10 



i; 









I 



x> 



z~5(age=1.2Gyr) 
z~7 (age=0.8 Gyr) 

— — Bouwens+07 
Bouwens+1 1 



J I I I IZj I I I I I lZi I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L. 



-23 



-22 



-21 
M 



-20 



-19 



1500 
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observations (blue filled circles) shows a moderate evolution at the bright end (Mt7v=-21.6 mag), consistent with previous observations. 
The open circles show our estimates if we include foreground contaminants namely brown-dwarfs. Thus, confusing brown-dwarfs for 
compact galaxies, which is possible for ground-based broad-band observations, will overestimate the UV LP. 



(Muv ~ -21.5). Thus, within -- 400 Myr the number 
density of Muv = —21 to —22 mag galaxies has increased 
by a factor of ~ 4 from z ~ 7 to z ~ 5, which imphes 
a rapid build-up of the star-formation rates in galaxies 
during this short period. 

4.3. Comparison to other z ^ 7 Samples 

While there are a few other studies focused on z ^ 7 
LBGs using HST observations, currently there are only 
four other ground-based surveys that have identified z ~ 
7 dropout galaxies using broad-band imaging (Figure 10: 
lOuchi et al.' '20091; iCastellano et al.l |2010'; 'Bowler et al.' 
1012; Hathi et al. I2012[ ). Using Subaru /Suprime-Cam 
in the Subaru Deep and GOODS-N fields, Ouchi et al. 
([2009. ) identified 22 z - 7 LBGs down to Y=26 mag. 
Our current survey reaches a similar depth with Ji=26.1 
mag (this filter is simil ar to the Y fi l ter), but covers a 
smaller survey area. Cast ellano et al.l ()2010[ ) found eight 
z-dropout galaxies in the GOODS-S field using HAWK- 
I and F0RS2 broad-band (Z, Y, J) observations in a 
survey area comparable to ours, how e ver w ith shallower 
survey depth. Recently iBowler et al.l()2012l) found about 
ten z ^ 7 LBG candidates in a large survey area (UltraV- 
ista) in the COSMOS field. While their survey depth is 
shallower (brighter than J=25 mag), their survey area 



is much larger allowing them to probe the brightest 
and ra rest galaxies at this redshift. Finally, Hathi et al.l 
(J2012D found two possible LBGs brighter than J=24.5 
mag from a slightly wider (~ 169 arcmin^) but much 
shallower (J - 25; 3a) survey in the GOODS-N field. 

On the other hand, using HST observations, sev- 
eral studies have found z ^ 7 dropout galaxies, all 
fainter (Mf/y> —21.5 mag) than the brightest galax- 
ies from ground-based observatio ns. All these studies 
use primarily t hree fields: HUDF (jBouwens et al.l 20111: 
IGrazian et al." '2011), GOODS-S (Wilkins et al." 'Ml 
iGrazian et al. 2011), or GOODS-N (Hathi et al. 2012|). 
While deep HST observations are appropriate in prob- 
ing the faint end of the UV luminosity function, ground- 
based surveys have the advantage of probing much larger 
volume and thus finding brighter and rarer galaxies, al- 
lowing us to constrain the bright end of the luminosity 
function as can be seen in Figure 9. All the observations 
brighter than Muv= -21.5 mag are currently obtained 
from ground-based observations. 

From Figure 9 we can see that at brighter mag- 
nitudes {MuY brighter than -22 mag), the best-fit 
Schechter fun ctions obtained for the HST observa- 
tions (Bouwen s et al.l [20071 120111 ) prefer a significant 
evolution from z ^ 5 to z ^ 7. On the other 
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hand, the fflound-b ased observations (|Ouchi et al.ll2QQ9t 
iBowler et alJl2Q12f ) suggest a less dramatic evolution at 
the bright end, however due to large error bars, strong 
evolution can not be ruled out. It is possible that in 
the absence of a clean method to identify contaminants, 
namely T-dwarfs (square symbols in Fig. 3), the ob- 
served z ^ 1 UV luminosity function that includes con- 
taminants (white filled circles in Fig. 9) will be over- 
estimate d and consistent with z ~ 5 luminosity func- 
tion (.Bouwens et al.|[2QQ7f ). This is especially likely to 
be the case for the ground-based broad-band imaging 
where atmospheric seeing makes it difficult to distin- 
guish between stars and unresolved galaxies based on 
their FWHM; zFourGE medium-band filters provide an 
effective method to identify T-dwarfs. Based on current 
observations, while the bright end of the UV luminos- 
ity function at z ~ 7 is reaching a better observational 
consensus, larger- area surveys with improved sensitivity 
are required to improve the UV luminosity function mea- 
surement and test for cosmic variance. 

5. PHYSICAL PROPERTIES OF Z - 7 GALAXIES 

Because near-IR medium-band filters provide a higher 
spectral resolution compared to the broad-band data, 
the constraints on physical parameters of these galaxies 
is improved (photometric redshifts, stellar masses, star- 
formation rates (SFRs), ages, extinction, etc.). Here, 
we constrain the physical properties of the two brightest 
galaxies, which have high S/N in the zFourGE bands. 

It is standard practice to compare galaxy photome- 
try with stellar population models that vary over some 
range of parameter space. Uncertainties on the stellar 
population parameters (e.g., mass, age, extinction, star- 
formation history, metallicity) are either derived through 
Monte Carlo methods or by marginalizing over some 
probability distribution function. Both of these meth- 
ods assume a prior that the models represent the data. 
In the sections that follow we compare the physical prop- 
erties derived using broad-band photometry alone, and 
then including the medium-band photometry. For illus- 
trative purpose, we plot the relevant figures either for 
the brightest (zF13965) or the 2nd brightest (zF25035) 
LBG (whichever yields extreme values) but tabulate the 
derived properties for the two brightest galaxies in Table 
2-3. 

5.1. Stellar Population Synthesis Modelling 

To derive the physical properties of our z ^ 1 galax- 
ies, we obtained model stella r population spectra using 
the iBruzual fc Charlotl (|2QQ3l ) stellar population synthe- 
sis code. We generated model spectra using a grid of 
dust, star formation histories (SFHs), and stellar popu- 
lation ages (t) ranging from 10 Myr to 750 Myr in loga- 
rithmic intervals. These limits on stellar age ensure that 
it is greater than the expected dynamical timescale of 
the galaxies but less than the age of the universe at a 
given redshift. We restrict the metallicities to Z=O.2Z0 
and assumed a Chabrier initial mass function (IMF) - 
for a Salpeter IMF the masses would be roughly 0.2 dex 
larger, but this would have no effect on the rest-frame 
UV-optical colors of the galaxies as the slopes of the IMFs 
at the high- mass end are identical between the two IMFs. 
In addition, we included nebular emission lines in our 
synthetic model spectra. These lines have been shown to 
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Fig. 10. — Comparison of photometric redshift probabihty dis- 
tribution of the brightest object (zF13965) using only broad-band 
photometry (i/S'T+Subaru+IRAC: hashed region), and adding 
medium-band filters (FourStar + i/S'T+Subaru +IRAC bands : 
filled red region). The photometric redshifts using the FourStar 
medium-band filters prefers a lower as well as narrower redshift 
range. 

affect the derived physical properties of galaxies in th e lo- 
cal universe (e.g., Papaderos et al. 2002; Pustilnik et al.l 
120041) and at high-redshifts (e.g., Schaerer & de Ba rrosI 
120091 : [Schaerer et al. 2011 : Fmkelstein et al. 2012, B. 
Salmon et al 2013 in preparation). We constrained 70 
to be negative (rising SFH), constant, or r >300 Myr 
as the SFHs of these galaxi es are expected (on a verage) 
to be increasing with time (Papovic h et al.l [2 01 If ). Next, 
we applied the dust attenuation t o our m odel spectra 
using the Calzetti dust law (Calzet ti et al.l 12000). which 
is appropriate to starburst galaxies. The spectra were 
then redshi fted and attenuated using the IGM attenu- 
ation from Meiksin (2006"). Finally, we computed the 
bandpass averaged fluxes in each of the filters. The best- 
fit model was obtained by minimizing the x^ between the 
bandpass averaged fluxes and the model spectra. To un- 
derstand the effect of different dust attenuation laws on 
the SED-derived parameters, we also used the SMC-like 
dust law (|Pei| [l992) described in Section 5.3. 

5.1.1. Photometric Redshift Distribution 

Figure 10 shows the photometric redshift prob- 
ability distribution for the brightest object. For 
the brightest object (zF13965), the photometric red- 
shifts derived for the iJST+Subaru bands tend to be 
higher with ^hot =8- 11^1.06' while inclusion of medium- 
band photometry lowers the photometric redshift to 
%hot=7.24lQ;25 (where the error bars are all 99% con- 
fidence; Figure 10). This is also true for the 2nd 
brightest object (zF25035), with Zphot=7.66lQ;57 and 

%hot=7.16lQ;^9 when derived using broad-band photom- 
etry alone, and then adding medium-band filters, respec- 
tively. The addition of medium-band photometry lowers 
the median redshift because of the Ji filter, which bet- 
ter constrains the Lyman-break. Surveys using a Y-band 

^'^ We model the star-formation rate as a function of e~*/'^. 
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Fig. 11. — Comparison of SED of our brightest candidate 
(zF13965) using all bands (top panel: if^ST+Subaru + FourStar 
medium-band +IRAC) and excluding medium-band photometry 
(middle panel: /fS'T+Subaru +IRAC bands). The SED model 
assumes O.2Z0 metallicity, a Chabrier IMF, the Bruzual &: 
Chariot 2003 model, and include nebular emission lines. The lower 
panel shows best-fit SED without IRAC bands. Adding FourStar 
medium-band filters improves the constraints on the photometric 
redshift and other physical properties, while exclusion of IRAC 
bands yields poorest constraints. 



near-IR filter may expect similar photometric redshift s. 
Including the medium-band filters, the photometric red- 
shift probability distribution functions favor lower red- 
shift solutions, with a narrower range of favored redsh- 
fits. We fit the stellar-population models to the measured 
photometry over a range of redshift spanned by these 
probability distribution functions, and we marginalize 
over these properties when determining constraints on 
stellar population parameters. 

5.1.2. Best-fit SED 

Figure 11 shows the best-fit model spectrum with and 
without the FourStar medium-band filters for our bright- 
est z ~ 7 LBG. The top panel shows the best-fit model 
spectrum using FourStar medium-band filters, HST and 
visible bands while the middle panel shows the best-fit 
SED without the FourStar medium-band data. The bot- 
tom panel shows the effects when the IRAC bands are 
excluded. In all three cases, we use the photometric red- 
shifts obtained from EAZY, using respective bands only. 
The behavior of the best-fit SEDs for the second bright- 
est z ^ 1 galaxy are very similar. 



5.1.3. UV Spectral Slope 

The UV continuum light is sensitive to dust absorption 
and its slope provides a good tracer of the dust attenu- 
ation for star-forming galaxies as the two are well corre- 
lated i n UV- luminous, star-forrning galaxies from z ^ 
to ^ 3 fMeurer et al."1995' ^T997',^T999';'Daddi et al."2QM 
Burgar ella et al. 2005; Laird et al. 2005; Reddv e t al.l 
2006 : [Dale et al.l 120071 : ITrever et al.l 120071 : ISalim etHI 

2001 . ^-^ i 

Recently, iFinkelstein et al.l (|2011[ ) measured UV spec- 
tral slopes oi z ^ 1 galaxies by directly fitting the best-fit 
SED with a power-law. We followed a similar procedure 
to measure /3 by fitting a power-law to the best-fit SED 
in the the rest-frame wavelength range 1250 < A < 2600 
A (|Calzetti et aLllTool . 

For the brightest LBG in our sample (zF13965), 
with this method we find /3=-1.88±0.1 using the HST 
bands while adding the medium-band photometry gives 
a steeper value of /3=-2.08±0.1 (see Table 3). For the 
2nd brightest LBG (zF25035), we find a similar trend in 
p with a shallower slope, /3 = -1.29±0.05 using HST 
broad-band photometry and a steeper P = — 1.44±0.12 
when adding medium-band photometry. For both galax- 
ies, our analyses favor steeper (i.e., bluer) UV spectral 
slopes when the medium-band photometry are included. 

The UV slope measured for both objects are between 
j3=-1.2 to /3=-2.1, within the range seen in other studies 
dX z r^ 7. Finkelstein et al (2011) found a median value 
of /3=-2.1 for galaxies with Muv = —20.6. Similarly, 
Dunlop et al (2012) found an average value of /3=-2 for 
bright z ^ 7 galaxies with Muv = —20.5. 

5.1.4. Constraints on Dust Attenuation from SED Fitting 

IFinkelstein et al.l (J2010D found the best-fit A(V)=0 - 
1.5 for a sample of LBGs at 2: ~ 7 obtain ed from the 
HST WFC3 in the HUDF field, while Labb e ^Fld] (l20Toh 
found A(V)=0 for a stacked sample of z ^ 7 LBGs, also 
obtained from the HST. Figure 12 shows the cumulative 
probability distribution of the dust content E(B-V) for 
our 2nd brightest object (zF25035) from the SED fit- 
ting. The shaded regions show 68% and 95% confidence 
regions. The top panel shows E(B-V) with broad-band 
+FourStar photometry, while the bottom panel shows 
the probability when including broad-band photometry 
alone (Table 3). 

While the best-fit model, with and without the 
medium-bands, yield similar range of E(B-V), the proba- 
bility distribution with the medium-bands spans a wider 
range with a median E(B-V) = 0.22, slightly lower 
compared with the median E(B-V) that excludes the 
medium-band filters. For the brightest LBG (zF13965), 
median E(B-V)=0.01 and E(B-V)=0.09 with and with- 
out FourStar bands respectively. For the reasons de- 
scribed above, the FourStar bands provide a more robust 
constraint on the physical properties of galaxies owing to 
their higher spectral resolution. 

5.1.5. Stellar Mass and Age 

Bowler et al. (2012) found that LBGs at ^ - 7 
with magnitudes in the range of our two most-luminous 
sources (J ^ 25 mag) have stellar mass ~ 5 x lO^M©. 
For fainter magnitudes, depending on the range of ex- 
tinction, the estimated stellar masses range from 3 x 10^ 
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TABLE 3 

Physical properties derived from spectral energy distribution fitting for the two brightest LBG candidates (zF13965 & 

ZF25035). 



ID 



Filters 



^phot 



/? 



E(B-Vi^ 



7 9/1+0.38 
'•^^-0.25 



+0.10 
0.10 
+0.10 
-0.10 
+0.05 
-0.05 



+0.06 
0.05 

+0.20 
0.01 



Log Stellar 
mass M0 



1^072- 

0.2 

+0.2 

0.3 

+0.5 

0.7 



Stellar age^ 
Myr 



Best-fit 
T (Gyr) 



ZF13965 if5'T+Subaru+IRAC+FourStar 
if5'T+Subaru+IRAC 
if5'T+Subaru++no IRAC+FourStar 



8.11 
7.24 



+0.45 

-1.06 

+0.38 

0.25 



-2.08 
-1. 

-1.17: 



0.09 
0.01 



9.8 
10.0 
11.3 



251 
159 
630 



100 
100 
100 



ZF25035 if5'T+Subaru+FourStar 7.16 

if5'T+Subaru 7.66 

if5'T+Subaru++no IRAC+FourStar 7.16 



^0735" 
0.19 

+0.77 
0.57 

+0.35 
0.19 



-1.44 
-1.29 
-1.52 



+07L2- 
0.12 

+0.05 
0.05 

+0.10 
0.10 



0.22 
0.24 



+T).cro" 

0.09 

+0.07 

0.20 



10.4 
10.2 
10.4 



+0T3" 
0.4 

+0.3 
0.3 

+0.5 
0.4 



321 
640 
640 



100 
100 
100 



' Median value. Errors indicate upper and lower 68% values. 

* We report the best-fit ages from the spectral-energy distribution modeling only, as the age constraints depend strongly on the assumed 
star-formation history. 
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Fig. 12. — Cumulative probability distribution of E(B-V), 
for the 2^^ brightest candidate (zF25035), using broad-band 
photometry alone (top panel), while the bottom panel includes 
medium-band photometry. While the median E(B-V)=0.2 are 
nearly same for both, the E(B-V) derived using medium-band 
filters allows a wider distribution compared to the the distribution 
obtained using broad-band photometry alone. 



- 1 X lQ^QM^(Finkels tein et alJl2QlQHLabbe"indl[2QTQl) . 
For our brightest object (zF13965), the favored range 
of stehar mass shifts to lower values after including the 
medium-band filters, with the most likely median mass 
decreasing by a factor of 1.5. As can be seen in Fig- 
ure 13, the probability distribution is narrower with me- 
dian mass 6.3 x 10^ M© when the medium-band filters 
are included; the absence of the FourStar medium-band 
photometry shifts the masses to higher values (median 
mass = 1 X 10^^ M©). On the other hand, for the second 
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Fig. 13. — Cumulative distribution of stellar mass for the 
most luminous z ~ 7 LBG (zF13965) with all bands included 
(top panel), excluding FourStar medium-band filters (middle 
panel), and excluding IRAC bands (lower panel). The addition 
of medium-band photometry shifts the stellar mass distribution 
to lower values (from 1 x 10^° M© to 6.3 x 10^ M©) because it 
better measures the spectrum of the galaxy. Removing the IRAC 
data allows for a large range of stellar-population parameters, 
including older age models, which greatly extends the range of 
stellar masses to higher values. 



brightest object (zF25035), the median mass is higher 
(2.5 X 10^^ M0) when medium-band filters are included. 
Figure 13 (bottom panel) shows the constraints if we 
remove the IRAC bands. In this case, the allowed range 
of stellar population parameters that fits the rest-frame 
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UV extends to models with very old ages and high mass- 
to-light ratios, and the fits are no longer able to exclude 
these possibilities. Thus, for the brightest object, the 
IRAC photometry is essential to better constrain the 
best-fit stellar masses. For the second object however, 
the presence of IRAC bands does not change the derived 
stellar mass (Table 3). A likely reason for this is that 
in our models the redshifted [OIII] line flux falls in the 
4.5/xm band over the photometric redshift range. This 
boosts the flux in the model, and requires lower stellar 
masses. 

The stellar ages for the two brightest objects range 
from 250 to 640 Myr, albeit with large uncertainties 
which depend on the assumed star-formation history. We 
can however estimate the timescale over which the stellar 
mass doubles. Because the best-fit r (Table 3) for both 
objects is large (100 Gyr), the timescale is simply stellar 
mass/SFR. This yields a timescale of about 150 Myr to 
double the stellar mass. 

In summary, using broad-band photometry alone, 
broadens the photometric redshift range compared to the 
result when the medium-band data are included. This is 
expected as the medium-band filters more accurately iso- 
late the location of the Lyman-break. In addition, cur- 
rently there are no broad-band filters corresponding to 
the central wavelengths of some of the FourStar medium- 
band filters (see Figure 1). 

While our comparison is based on only two z r^ 7 
dropout galaxies, they are at the upper luminosity (and 
stellar mass) range of galaxies at this redshift. For our 
small sample, adding the medium-band photometry to 
the SED fitting yield more accurate constraints on the 
physical properties of the z ^ 7 galaxies compared to 
using broad-band photometry alone. To generalize these 
conclusions to the z ~ 7 galaxy population will however 
require larger samples with deep medium-band imaging. 

5.1.6. Effect of Lyman-a Line and Ji Filter on SED 

Fitting 

We further investigated the presence of a Lyman-a 
line and the effect of the Ji filter on the derived phys- 
ical properties of our bright LBGs. To do this, we ob- 
tained the best-fit SED (section 5.1) by first including no 
Lyman-Qf emission and then by excluding the Ji filter. 
Excluding the Ji filter or absence of a Lyman-a flux 
results in essentially the same values as that for broad- 
band only derived parameters. This reinforces the idea 
that we need Lyman-a equivalent width measurements 
to understand (1) the dust extinction and (2) the stellar 
masses, especially for galaxies with Lyman-a emission 
line. 

5.1.7. Lyman-a Equivalent Width 

In principle, the Lyman-a equivalent width W pro- 
vides a proxy to the star-formation rate in the galaxy, 
and defined as W^ = F\ine/F^^ where F^ne is the total 
line fiux (erg s~^ cm~^) and F^ is the continuum flux 
density (erg s~^ cm~^ A~^). If our z ^ 7 LBG candi- 
dates have Lyman-a line, it will contribute a flux excess 
in the Ji filter, the reby decreasing (i.e., m aking brighter) 
the Ji magnitude (|Papovich et al I [200 If ): 
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TABLE 4 

Estimated rest-frame Lyman-q; equivalent widths and 

LyMAN-Q; line fluxes FOR THE TWO BRIGHTEST Z r^ 7 LBG 
CANDIDATES. 



ID 



EWo Lyman-a flux 

(A) (IQ-^^ergs-^cm-') 



Method 1 

ZF13965 

ZF25035 



16.5^ 



:; + 15.3 
-15.3 

26.2ll^e:l 



1.15Z 



' + 1.07 
-^-1.07 
1 07+0.66 
-l--U< _0.66 



Method 2 

ZF13965 

ZF25035 



-1.0^ 



^+12.6 
-12.6 

10 c: + 14.3 
-LO-0_14.3 



0.0: 



-1+1.0 
-1.0 

^•'-'-0.66 



where Wq is the rest-frame Lyman-a equivalent width 
and AA is the Ji filter width. 

To measure Atti, we followed two different methods to 
take advantage of the presence of Ji filter. In the 1^* 
method, we first obtain the best-fit SED using all pho- 
tometry and including the emission line (see Section 5.1). 
We then use this best-fit SED, but remove the emission 
hue, and measure the bandpass Ji magnitude. This mag- 
nitude compared with the observed Ji magnitude yields 
Am. In the 2"^ method, instead of removing the emis- 
sion line, we obtain the best-fit excluding the Ji filter and 
compare the observed Ji magnitude with the synthesized 
Ji magnitude to get Am. We calculate the rest-frame 
equivalent widths using Equation 8 for our two brightest 
z ^ 7 LBG candidates (Table 4). Using these Wo, we 
calculate the Lyman-a line flux: 



Fune = F^Woil+z). 



(9) 



(8) 



Table 4 shows the estimated Lyman-a EWs and 
line fluxes, obtained from both methods, for the the two 
brightest z ^ 7 LBG candidates. The Lyman-a EWs 
are smaller when the Ji filter is excluded ( 2"^ method). 
On the other hand, the Lyman-a fluxes obtained from 
the 1®* method are similar to the recent observations 
of spectroscopically confirmed LBGs at 2^=7.008 and 
z =7.109, with line fluxes 1 62 xlQ"^^ and 1.21 xlO"^^ 
erg s~^ cm~^, respectively (Vanzell a et al.|[2QTll ). These 
line fluxes are accessible to deep spectroscopy with >8m 
telescopes, and thus these two candidates provide good 
targets for future spectroscopic followup. 

5.2. Specific Star Formation Rate 

The specific star formation rate, sSFR (SFR/M^) pro- 
vides a way to quantify the evolution of mass build-up 
given a certain SFR. Many observations suggest that the 
sSFR is nearly cons tant from z ^ 7 to z ^ 4 (Stark et al.l 
[2QQ9I : iGonzalez et al.l l2QlQl : iBouwens et all 12011). This 
implies that the SFR increase from high-redshifts to 
low-redshifts, in conflict with the general assumption 
of constant or declini ng SFH. However, rece nt studies 
(e.g.. Star k et al.l[2012l : ISchaerer fc de Barrosll2012f ) sug- 
gest that there is a stronger evolution in the sSFR from 
low-redshift to high-redshifts, than suggested by previous 
studies. Thus, to resolve this uncertainty among differ- 
ent studies, large samples of high-redshift galaxies are 
needed. 

Here we estimate the sSFR for the two brightest galax- 
ies in our z ^ 7 sample using the more accurate physical 



16 



Tilvi et al. 



10,0 r 



o 



q: 



(/) 
en 



KOr 



OJ r 









1 I 



I 

o 



A 



* 



I 



A 



©This work 
AGonzalez+10 



OSalim + 07 

■ Noeske+07 

■ Daddi+07 
^Reddy+12 

■ Bouwens+1 1 
O Stark+09 



A(1500A) [mag] 
2.0 4.0 6.0 



2 4 6 8 

Redshift 

Fig. 14. — The redshift evolution of the specific SFR. Blue 
filled circles represent the two brightest galaxies from FourStar 
observations. Open diamond, open circles, and open triangle 
represent data from Salim et al. (2007), Stark et al. (2009), and 
[Gonzalez et al. (|2010i ) respectively while filled symbols- grey 
squares, cyan s quar e, sta rs, and red squares show observations 
from Noeske et aP (12007.). .Daddi et al. (2002,), .Reddv et all 
([2012), and Bouwe ns et al.l ((2^011) respectively. These galaxies 
have stellar masses M^ = (1 — 2) x lO-*^*-* Mq. The error bars on 
the FourStar observations are the uncertainties in the estimated 
stellar mass obtained from the best-fit SED. The sSFR remains 
nearly constant over 1 Gyr period from z = 8 to z = 4, which 
favors the i dea of rising star-for mation history with increasing 
stellar mass (jPapovich et al.ir2011l) . 

parameters derived from the fits to the FourStar medium- 
band + HST + IRAC data. We derived the SFRs us- 
ing the best-fit SEDs: for the brightest LBG (zF13965), 
the median SFR=77.6l42 M© yr~^, much smaher com- 
pared with the 2nd brightest LBG (zF25035) with me- 
dian SFR=285.7l?^^ Mq yr-^ (Table 5). 

Our derived sSFR '^13 Gyr~^ for the two brightest 
z ^ 7 LBG candidate is shghtly higher compared with 
the values derived for lower luminosity (and lower mass) 
galaxies at z ~ 7, which implies a nearly co nstant sSFR 
from z ^ 7 to z ^ 4 (e.g., Reddy et al.l 1201 2) for galaxies 
over approximately a decade in stellar mass (2 x 10^ Mq 
to 2 X 10^° Mo). These sSFRs favors the idea that galax- 
ies at these epochs have rising star-formation histories 
such that the sSFR is nearly independent of mass at a 
fixed redshift (e.g. Papovich et al. 2011). Models with 
constant or declining SFHs would imply decreasing sSFR 
with increasing mass at fixed redshift, which is disfa- 
vored. 

5.3. Comparison Using Calzetti & SMC Dust 
Extinction Laws 

To compare the effect of different extinction laws on the 
derived physical properties of 2; ~ 7 LBGs, we derived the 
best-fit SEDs (see Section 5.1.2) and physical properties 
of the two brightest z ^ 7 LBG candidates using Calzetti 
and SMC dust extinction laws. 

For the brightest LBG (zF13965), all derived param- 
eters (stellar mass, age, and x^) using SMC dust ex- 
tinction law decrease except for the dust extinction (see 
Table 5). For the 2nd brightest LBG (zF25035), stel- 
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E(B-V) 



0.6 



Fig. 15. — Comparison of dust extinction derived from the best- 
fit SED using Calzetti and SMC dust extinction law. Here we have 
shown the dust extinction of the 2nd brightest LBG (zF25035) as 
this object yields significantly different extinctions for SMC and 
Calzetti dust extinction law, compared with the brightest LBG 
(zF13965). The SMC-like law yields a median E(B-V)=0.1, a value 
much smaller compared with median E(B-V)=0.22 obtained using 
the Calzetti dust extinction law. 



TABLE 5 
Comparison of derived physical parameters from the 

BEST-FIT SEDs USING CALZETTI AND SMC DUST EXTINCTION 
LAWS FOR THE TWO BRIGHTEST LBGS. 





ZF13965 


ZF25035 


xV^- 


2.3 [1.8] 


0.7 [0.7] 


Log Mass(M0) 


9-8tS:2 [9-6t°o:4] 


io.3tS:^, [9.8tS:^j 


E(B-V) 


omtVol [0.05lS:°o3] 


p. r)r) + 0.09 Ip. -1+0.05] 


SFR(M0yr-^) 


77.6tll:l [57.8lg:I] 


285.7lf/,:« [62.7lJ|:^] 


sSFR(Gyr-^) 


13.5 tl\l [13.5 tl\j] 


13.5 tr, [10.7 tli] 



Values in square brackets are derived using SMC extinction law. 



lar mass, extinction and x^ either decrease or remain 
unchanged except for the stellar age. A major differ- 
ence among these two LBGs is that the dust extinction 
for zF25035 is significantly reduced (Auv=0.9) when us- 
ing the SMC extinction law (see Fig. 15). This smaller 
extinction value also gives rise to the lower SFR com- 
pared with the Calzetti dust extinction-derived SFR. 
This demonstrates that the assumption of the extinc- 
tion law affects the interpretation of the dust content in 
high-redshift galaxies. In particular, the SMC-type law 
- which may be mor e applicable in these distant galaxies 
(|Oesch et al.ll2Q12bl ) - reduces the implied dust content 
compared to a Calzetti- type law for low-redshift lumi- 
nous starburst galaxies. 
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6. SUMMARY & CONCLUSIONS 

We have obtained a sample of z ^ 7 galaxies from 
a deep near-IR survey (zFourGE ) using medium-band 
filters for the first time. We define a color-color selec- 
tion criteria for the near-IR medium-band filters, which 
cleanly isolate z ^ 7 LBG candidates from brown dwarf 
stars and low-redshift galaxies. Using our criteria we 
have identified three robust z ^ 7 candidate LBGs in a 
survey area of about 155 arcmin^ in the COSMOS field, 
which has very deep publically available data. The avail- 
ability of near-IR medium-band and broad-band photom- 
etry allowed us to compare the derived physical proper- 
ties of high-redshift galaxies and how these properties are 
influenced by the choice of filter widths. Our principal 
findings are summarized below. 

Stars vs compact galaxies— While the contamina- 
tion from nearby T-dwarf stars is a serious concern 
for dropout-selected galaxies, especially when using 
the ground-based broad-band photometry, the FourStar 
medium-band filters allow us to cleanly distinguish be- 
tween nearby dwarf stars and high-redshift galaxies that 
are relatively compact. 

One of the three z ^ 7 LBGs in our sample (zF25035) 
is very compact and indistinguishable from a point source 
at ground-based resolution (and possibly unresolved even 
by HST). Our ability to distinguish this object from cool 
stars is possible only due to the availability of multi- 
wavelength medium-band photometry. Moreover, this 
demonstrates that there exists a population of very com- 
pact galaxies at z ~ 7. 

Bright end of the z ^ 7 UV luminosty function— Us- 
ing the number density of z ^ 7 galaxies obtained us- 
ing dropout technique and corrected for incompleteness, 
the rest-frame stepwise UV luminosity function at z ~ 7 
shows a moderate evolution from z ^ 5. The number 
density of bright LBGs (M[/y~-21.5) increases by a fac- 
tor of 4 from z ^ 7 to z ^ 5 albeit with large error 
bars. This is consistent with other ground-based studies 
but implies lesser evolution compared to the factor of 10 
increase implied by some iJST-based studies. 

Physical properties of z ^ 7 galaxies— In general, the 
presence of medium-band photometry yield tighter con- 
straints on the photometric redshifts and improved phys- 
ical constraints on the stellar population parameters of 
these galaxies. 

(i) For our two brightest LBGs, the best-fit SED 
model with broad-band + medium-band photom- 
etry yields lower photometric redshifts, and much 
tighter redshift probability distribution compared 
to the %hot~8 derived using only broad-band pho- 
tometry (Figure 10). This is also true for our 2nd 
brightest object. 

(ii) For both objects, the UV spectral slope P tends to 
be shallower when derived using broad-band pho- 
tometry alone. For the brightest object, P changes 
from -1.88 to -2.08 when FourStar medium-band 
photometry is included. 

(ii) The SED modehng prefers a lower dust attenuation 
with median E(B-V)=0.01 for broad-band photom- 
etry, while it increases to median E(B-V)=0.09 
when including medium-band photometry. The 
FourStar bands allow for redder UV colors which 



increases the range of permitted values for E(B- 
V). Comparing the above values, derived using the 
Calzetti dust law, with the E(B-V) derived using 
SMC-like dust law, we find that the probability dis- 
tribution from the latter case is much narrower and 
shifts towards lower values (Figure 15). 

(iii) Including the FourStar medium-band filters, the 
SED modeling yields narrower probability distri- 
bution for stellar mass for the two most luminous 
objects in our sample. We argue that adding ad- 
ditional bands shifts the allowed regions of the 
stellar-population parameter space, improving our 
physical understanding of the galaxies. We find 
that the presence of IRAC photometry is also 
important to better constrain the best-fit stellar 
masses. 

(iv) We predict the Lyman-a equivalent widths and 
line fluxes of the two brightest candidates and find 
that the presence of Lyman-a line influences the 
derived physical properties of z ^ 7 candidates. 
The estimated rest-frame equivalent widths (17 and 
26 A) make these two objects good targets for fu- 
ture spectroscopic followup. 

(v) The sSFRs for the two most luminous galaxies 
are slightly larger compared with other lower mass 
galaxies at z ~ 7. This favors the ideal that the 
SFR increases with increasing stellar mass at this 
redshift. 

While it is early to reach robust conclusions about the 
effect of medium-band photometry on the derived prop- 
erties of z ^ 7 galaxies, medium-band filters will likely 
provide better constraints on the physical properties, due 
to their narrower redshift probability distributions, com- 
pared to the broad-band photometry. If this is indeed 
true, in the absence of spectroscopic redshifts, multiwave- 
length medium-band photometry will provide the best 
constraints on the physical properties of z ^ 7 galaxies. 
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